hep-ph/0009003 



(N 

O 
O 
(N 

< 



9" 

43 



ctf 



A+ Production in Polarized pp Scattering and 
Polarized Gluon Distribution 



Kazumasa OHKUMA 1 )'^, Kazutaka SUDOH ^'Q and Toshiyuki MORII 1)i2), Q 



1) Graduate School of Science and Technology, Kobe University 
Nada, Kobe 657-8501, JAPAN 
(N, 

2) Faculty of Human Development, Kobe University 

^ i Nada, Kobe 657-8501, JAPAN 

O 

o 
o 
o 



Ph! abstract 

To extract information about the polarized gluon distribution, AG(x, Q 2 ), in 
• i— i , ^ 

the nucleons, we propose A+ productions in polarized pp scattering, p+p — > A++X, 
H ' 

which will be observed at RHIC experiment starting soon. For this process, we 

have calculated the spin correlation differential cross section, dAa/dp T , and the 

spin correlation asymmetry defined by A LL = [dAcr/dp T ]/[da/dp T }. We have found 

that the All is sensitive to the polarized gluon distribution in the nucleon and thus 

the process is promising for testing AG(x, Q 2 ). 
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The proton spin has long been considered to be given by the sum of the 
constituent quark spin, AS, in the naive quark model, which suggests AS = 
Am + Ad = 1 and As = 0, where Am, Ad and As denote the amount of the proton 
spin carried by u, d and s quark, respectively. However, the experimental data on 
polarized structure function of nucleons, g{ (x, Q 2 ), have revealed that the con- 
tribution of the constituent quark spin to the proton spin is quite small, AS ~ 0.3, 
and the s quark is negatively polarized with a rather large value As ~ —0.12 Q. 
This has been called the "proton spin puzzle" ||. Actually, the proton spin is 
given by the sum of the spin of the quarks (valence and sea), gluons and the 
orbital angular momenta among them; 

i = ^AS + AG + (L z ) q+9 , (1) 

where AG and (L z )q+ g represent the amount of the proton spin carried by the 
gluon and the orbital angular momenta of quarks and gluons, respectively. To 
understand the physical ground of the sum rule of Eq. ([!]), it is important to know 
the behavior of the polarized partons in the nucleon. However, knowledge about 
the polarized gluon density, AG(x, Q 2 ), in a nucleon is still poor, though many 
processes have been proposed so far to extract information about it. On the other 
hand, the RHIC experiment which will start soon, will provide us many useful data 
for extracting information about polarized gluons in the near future. 

In this paper, to extract the polarized gluon distribution AG(x,Q 2 ), we pro- 
pose another process, p + p — > A+ + X (Fig. [I]), which would be observed in the 
forthcoming RHIC experiment. In this process, A+ is dominantly produced via 
fragmentation of a charm quark originated from gluon-gluon fusion. The reason 
we focus on this process is as follows. The A+ is composed of a heavy quark c and 
antisymmetrically combined light u and d quarks. Hence, the A+ spin is basically 
carried by a charm quark which is produced via gluon-gluon fusion at the lowest 
order in this process as shown in Fig. |2| Q. Therefore, observation of the spin of 
the produced A+ gives us information about the polarized gluons in the proton. 
In order to study this attractive process; 

"- 1 Since charm quarks are not main constituents of the proton, the gluon-gluon fusion process 
is dominant for charm quark production. 
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Figure 1: The diagram for p + p 
A+ + X at the lowest order. 



Figure 2: The subprocess at the low- 
est order. 



p(Pa)+P(Pb)-+K(Pa+)+X, 



(2) 



whose subprocess at the lowest order is 



g(p a ) + giPb) -» c(p c ) + c(p 5 



(3) 



where p t denotes the four-momenta of the i-particle " 2 and the over-arrow means 
that the particle polarization is either provided (g) or measured (c), we introduce 
two useful observables which will be measured by the forthcoming RHIC experi- 
ment: one is the spin correlation differential cross section, dAa/dp T , and the other 
is the spin correlation asymmetry, All- They are defined as follows; 

dAa _ dcr(++) - dcr(-\ — ) + dcr( ) - da( — h) 



dp T 
A LL = 



dp T 

[da(++) - da(A — ) + da(- 



(4) 



da( — h)] /dp q 



[da(++) + da(+-) + da( — ) + da(-+)] /dp 7 
dAcr/dp T 



da/dpj 



(5) 



where da(-\ — )/dp T , for example, denotes the spin-dependent differential cross sec- 
tion with the positive helicitiy of the target proton and the negative helicity of the 
produced A+. 

Let us consider the process in the proton-proton cm. frame and take the 
four-momenta of p % as follows; 



Pa 



^-(l,±/3,0) with /3 = yi-_ £ 



" 2 For example, if i = A+, pi indicates the four-momenta of A^ 
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Pa+ = ( e a+,Pl>Pi 



'm 2 A+ + p 2 cosec 2 Q,p T cot Q,p T ), 

Pa, b = X a. b P A ,B , Pc = —1 (6) 



where the first, second and third components in parentheses are the energy, the 
longitudinal momentum and the transverse momentum, respectively 0. and m* 
represent the cm. scattering angle of the produced A+ and mass of the i-particle, 
respectively. x a ^ and z are the momentum fraction of the proton carried by the 
gluon and the one of the charm quark carried by A+, respectively. Notice that 
we assume that the scattering angle of the A+ produced in the final state and the 
angle of the charm quark in the subprocess are almost same, 8 C ~ 8 A + = 0. This 
assumption is not unreasonable, because the momentum of A+ is almost carried 
by a charm quark whose mass is much larger than the other constituents of A+. 
Then, using Eq.(|6|), we define the following Lorentz invariant variables; 

s = s - 2m 2 , 



t = t — m — m.+ = — y/s \/m 2 . + + p 2 cosec 2 9 + j3p cot 



p ~ m A+ = ~^ \/ m A+ + Plcosec 2 Q - (3p T cot , (7) 



■■~aj v ~ |_y ■■- A + ■ r T - 

u = u — mi ' /_ 

where s, t and u are conventional Mandelstam variables, s = (p A + p B ) 2 , t = 
(p A + —p B ) 2 and u = (p A + — p A ) 2 , respectively. Furthermore, for the subprocess, 
we define 

t\ =t — m c = — t, U\ = u — m c = — u, (8) 

with the subprocess Mandelstam variables s = (p a + p b ) 2 , t = (p c — p h ) 2 and 

Using Eq.(|^)~Eq.(Rp, the spin correlation differential cross section(Eq.([|)) can 
be expressed as 

dAa /-© max r 1 r 1 



IT = L. jL L a G PA ^ 9a (x a ,Q 2 )AG fB ^ b ( Xh Q 2 )AD s ^(z) 

"T a b 

x — —Jdx a d.XbdQ, (9) 

dt 



" 3 Hcrc we do not neglect the masses of all particles, though we basically follows the method of 
Rcf. 1. 
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where G PA ^ ga (x a , Q 2 ), AG^ B ^ b (xb,Q 2 ) and AD^+(z) represent the unpolar- 
ized gluon distribution function, the polarized gluon distribution function and the 
spin-dependent fragmentation function of the outgoing charm quark decaying into 
a polarized A+, respectively. Moreover J is the Jacobian which transform the vari- 
ables z and t into and p T . Unfortunately, at present there is no established spin- 
dependent fragmentation functions because of lack of experimental data. However, 
since a charm quark is much heavier than other constituents of A+, it is expected to 
be very rare for a charm quark to change its spin alignment during fragmentation 
process. Therefore, it is not unreasonable to substitute D C ^ A + for AD^^+ . In this 
work, we use the model by Peterson et al. (3| for both D C ^ A + and ADg_^+ F° r 
the subprocess, the spin correlation differential cross section, dAa/dt is calculated 
to be 



dAa ira 2 



dt s 

and the Jacobian, J, is given by 



m 2 j 9ti - 19ui 8s 1 s J ti - Mi 1 3 f 2t x 



(10) 



where z is 



2s/3p 2 cosec 2 . , 

j= _ J2L T (li) 

zs x m 2 + + p 2 cosec 2 



z = ^ + ^ (12) 

X a Xf) 

with X\ = —t/s and X2 = —u/s. The minimum of x a , Xb are given by 

min X\ mill X a X2 , > 

x = , x. = . (13) 

1 — X 2 X a — X\ 

The unpolarized differential cross sections were calculated by Babcock et al. ||. 

For numerical calculation, we use as input parameters, m c = 1.5 GeV, m p = 
0.938 GeV and m A + = 2.28 GeV |J. We limit the integration region of and 
p T of produced A+ as § < < $* and 3 GeV < p T < 20 GeV for ^ = 200 
GeV and 500 GeV, in order to get rid of the contribution of the diffractive A+ 
production and also the A+ production through a single charm quark production 
via W boson exchange and W boson production. As for the gluon distributions, 



" 4 Actually, we used the Peterson functional form presented by the Particle Data Group 
which is normalized as jD c _ rA +(z)dz = 0.503. 
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we take the GS96(set-A and -B) [[?]] and the GRSV96 || parameterization mod- 
els for the polarized gluon distribution function and the GRV95 || model for the 
unpolarized one. Note that the GS96 and GRSV96 models can excellently repro- 
duce experimental results for the polarized structure function of nucleons, though 
behavior of these polarized gluon distributions is quite different. In other words, 
the data on polarized structure functions of nucleons and deuteron alone are not 
enough to distinguish the model of gluon distribution functions. We are interested 
in the sensitivity of those observables of Eq. (§) and Eq. (|5|) on the polarized gluon 
distributions in the nucleon in this process. 

We show the p T distribution of dAo/dp T and An in Fig. |3] for y/s = 200 GeV 
and in Fig. [| for y/s = 500 GeV, respectively. Notice that in Fig. |] the absolute 
value of dAo jdp T is presented, since the negative value of dAo jdp T cannot be 
depicted in the figure which has an ordinate with logarithmic scale. Actually, 
for the case of y/s = 500 GeV, the value of dAa/dp T becomes negative for p T 
smaller than the value corresponding to the sharp dip shown in Fig. |j. On the 
other hand, dAo/dp T at y/s = 200 GeV is positive for all p T regions as shown in 
Fig. |3|. To understand this quite different behavior of dAo jdp T depending on y/s, 
some comments are in order. As seen from Eq.(^)~ Eq. flTB|) , the sign of dAo /dp T 
is determined by the sign of dAo/dt because all variables except for dAo/dt are 
positive in whole kinematical regions of x a , Xb and O. Since both t\ and u\ are 
negative and smaller than s in magnitude, the sign of dAo/dt mainly originates 
from the second term of the right hand side of Eq . (|10D , in which t\ —U\ is calculated 
to be 

i x — «i ~ x b y/s I X a yf~S \\J m \+ + PyCosec 2 + p T cot Qj [ , (14) 

in the limit of the massless proton 0. Furthermore, in the same limit, x™ in reduces 

to 

\/m 2 . + + p 2 cosec 2 6 + p T cot 

c * — v , ; — x- (is) 

y/s — ( ,/m 2 + +p 2 cosec 2 — p T cot 



" 5 In the region of s much larger than nil where we are now focusing, this limit with s = s and 
j3 = 1 is a good approximation. 
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Figure 3: The unpolarized and spin correlation differential cross section (left panel) 
and the spin correlation asymmetry (right panel) as a function of p T at y/s = 200 
GeV. The solid line in the left figure represents the unpolarized differential cross 
section with the GRV95 model for the unpolarized gluon distribution. The dashed, 
long-dashed and dot-dashed lines indicate numerical results with the set A, B of 
the GS96 model and the GRSV96 model, respectively, for the polarized gluon 
distribution. The error bars for the dashed line shows the statistical sensitivity 
(see text). 
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Figure 4: The same as in Fig. y, but for \fs = 500 GeV. 



from Eq . flilf ) . Eq. flTBl) shows that x a becomes very small in some region of p T and 
O when ■y/s is large and hence leads to large values of the gluon distributions, 
G PA ^ ga (x a ,Q 2 ). Thus when Eq . (|14"D becomes negative due to very small x a , the 
negative contribution of dAa/dt to dAa/dp T becomes larger for ^/s = 500 GeV 
than for y/s = 200 GeV because G PA ^ ga (x a , Q 2 ) is larger at smaller x a . Therefore, 
for large i/s such as y/s = 500 GeV, even after integration, dAa/dp T is negative 
for p T smaller than the value corresponding to the dip, as shown in Fig. ^. It is 
very interesting to note that the behavior of A LL strongly depends on y/s as shown 
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from Figs. |3| and [|. In any case, both figures show that All is sensitive to the 
polarized gluon distribution function. 

Finally, to examine if our predictions can be tested at the forthcoming RHIC 
experiment, we estimate the statistical sensitivity of the spin correlation asymme- 
try, 5 All, according to the method given in Ref (TIJ . The value of 5 All is estimated 
by 

1 l 

P \/ b A+ £ £T a 



We can estimate the statistical sensitivity, 5 All, for T = 100-day experiments by 
using the parameters of the beam polarization (P = 70%), a luminosity (£ = 
8 x 10 31 (2 x 10 32 ) cm^sec" 1 for y/s = 200 (500) GeV), the trigger efficiency 
(s = 10% ) of detecting high p T charm production events through their semi- 
leptonic decays and a branching ratio (b A + = _Br(A+ — ► pK~n + ) ~ 5% ffl) Q a 
denotes the unpolarized cross section integrated over suitable p T region. The results 
are summarized in the right panels of Figs. |3] and ^. As shown in these figures, 
5 All is smaller than the difference of the model predictions at moderate p T region 
and thus our predictions are expected to be tested in the RHIC experiment, though 
the differential cross section, da/dp T , becomes rapidly smaller with increasing p T , 
and thus in the larger p T region, 5 All becomes too large to distinguish the model 
of polarized gluon distribution functions. Furthermore, note that 5 All is expected 
to become small if a trigger efficiency e is improved and/or another decay modes 
of A+ are additionally taken into account. Moreover, the statistics should be twice 
by using the A^ data as well, since high p T charm production is dominated by 
gluon fusion, as long as the A+ decays satisfy CP invariance. 

In summary, to extract the polarized gluon distribution AG(x, Q 2 ), we pro- 
posed a new polarized process, p + p — > A+ + X, which could be observed in the 
forthcoming RHIC experiment. We calculated the the spin correlation differential 
cross section, dAa/dp T , and the spin correlation asymmetry, All, and found that 
All is quite sensitive to the polarized gluon distribution functions; we can rather 

" 6 Here, we simply assumed the efficiency for reconstructing the pK~ir + decay of A+ to be 100 
% by quoting only the decay branching fraction. Although the actual efficiency might be more 
likely to be less than 10 - 20 % rather than more than 50 % (almost perfect, background free 
measurement), this should be studied by experimentalists in the forthcoming RHIC experiment. 
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clearly distinguish the model of polarized gluon distribution functions at moderate 
p T region. It is also remarkable that the All shows very different behavior in the 
region of ^/s = 200 ~ 500 GeV, covered by RHIC Therefore, the process looks 
promising for testing the model of polarized gluon distribution functions, though 
the present analysis is based on the leading order calculation. In order to get more 
profound information on the behavior of polarized gluons, we need the next-to- 
leading order calculation. Furthermore, we also need further investigation of the 
polarized fragmentation functions, AD S _^^+ (z) , to get more reliable prediction. Al- 
though these subjects are interesting and important in their own right, they are 
out of scope in this work. 

We hope our prediction will be tested in the forthcoming RHIC experiment. 
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